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ABSTRACT 

Continued studies were made of bluff-body separation and its 
effect upon the response of an aeroelastic model for a range of tunnel 
air speeds and model cross sections. The motivation for this work 
was the application to the antenna- element vibration problem as posed 
in the various Wullenwebber type antenna installations operated by the 
United States Navy throughout the world. 

A retrofit configuration change was devised to significantly 
decrease the wind-induced response of the long horizontal '‘^boom boards". 
Although flow separation stS.Il exists, the character of the aerodynamic 
forcing function has been modified such that the energy spectra have 
apparently been broadened from a single dominant peak into many weaker 
peaks. Thus the effective total response is reduced. 

The digital data processing has been extended in scope due to the 
repeatability in signal analysis. Estimates have been made of the "two- 
dimensional" random input forcing function using the output response in 
conjunction with the aeroelast ic-model transfer function. Reasonably 
good results were obtained and an estimate was made of the root-mean- 
square value of the unsteady fluctuating air load. 
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1. Introduction. 

A bluff body may be defined from an aerodynamic sense as a body whose 
drag forces are primarily due to pressure or form drag. The airflow does 
not progress smoothly around the body, but instead separates randomly and 
defines a turbulent flow region. As a consequence of the flow separation, 
fluctuating pressures or air loads are generated which can be resolved 
into unsteady air load components in both the lateral (lift) and stream- 
wise (drag) directions. In general, the fluctuating lift forces dominate 
and give rise to the most significant vibration modes. 

Physical structures are elastic and consequently will respond to 
dynamic air loads that are wind generated. Examples of the response of 
bluff bodies to unsteady air loads are frequently observed in every day 
events, examples of which would include the galloping transmission line 
problem, vibrating telephone wires, smokestack motion, ground-wind loads 
on missiles, and vibrating -antenna elements. There are many other 
situations of structural shapes responding to wind-induced air loads. 

It was the purpose of this study to investigate the bluff-body 
problem with particular application to an important element of a large 
antenna (Wullenwebberttype) system that is currently in use by the 
United States Navy. The basic arrangement includes a number of towers 
located in a circular pattern. Figure 1 shows a pair of towers supporting 
a horizontal boom board. The boom boards extend around the entire circle, 
and serve the function of supporting vertical reflecting wires. The boom 
boards are 28 to 56 feet long (depending on the installation) and are 
approximately 100 feet above the ground. There have been instances 
reported of large amplitude motion of boom boards which were directly 
exposed to the airflow, and even some instances of component failures. 

An important consideration in this study is that the cross sections 
used for the boom boards are not the same for every installation since it 
is the option of each contractor to perform the detail design. Figure 2 
is a sketch showing the basic cross sections used for the boom boards. 

Only one installation (Winter Harbor, Maine) employs the circular cross 
section. There are three basic types of rectangular cross sections. The 
solid laminated boom board is used in three general sizes. The hollow 
square and the inverted *'U" are lighter per unit cross section area, but 
lack the strength of solid construction. The inverted does have 
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short lengths of wood across the bottom which support the vertical wires 
and provide structural support to maintain the shape. 

In order to fulfill the objectives of this study, it was desirable 

to : 

(1) Gain an insight into the factors influencing the response of a 
simulated aeroelastic bluff body. 

(2) Determine possible fixes suitable for application to the actual 
antenna installation. 

(3) Attempt to estimate the character and magnitude of the aerodynamic 
input . 

This work has been sponsored by the Structures Division of the 
United States Naval Civil Engineering Laboratory (USNCEL), Port Hueneme , 
California under work request WR 600028, dated 23 August 1965. 

The investigation was conducted with Lieutenant R. H. Jesberg, 

U. S. N. Tests were conducted in the West Coast Research Corporation 
(W.R.C.) low-speed wind tunnel located in Building 234 at the United 
States Naval Postgraduate School, Monterey, California. For detailed 
information on equipment and data handling procedures used in this study, 
see "Data Handling Procedures in Support of the Study of Aeroelastic 
Problems in an Antenna Vibration Project" by Lieutenant R. H. Jesberg, 

U. S. N. [7J. 

Sincere appreciation is extended to faculty and staff of the 
Department of Aeronautics for assistance and support which was vital to 
the success of this investigation. The author is indebted to Associate 
Professor L. V. Schmidt for his invaluable counseling and guidance. 
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2. 



Test Description. 

An aeroelastic bluff body model was installed to vertically span 
the 3.5 X 5.0-foot test section of the W.R.C. low-speed wind tunnel. The 
wind tunnel is of the continuous flow type, and has an atmospheric vented 
test section. Test parameters included tunnel dynamic pressure (q), 
model cross-section geometry, and model natural frequency 

The cross section of the model was altered by fitting various shapes 
of balsa wood on an aluminum I beam. The low density of balsa wood made 
it possible to make changes in cross-section geometry without signif- 
icantly altering the weight, center of gravity, and mass moment of inertia 
of the. complete model. 

The aluminum I beam was flexurally supported at the top (outside of 
the test section) and could move laterally in response to the unsteady- 
lift loads. A drag constraint was provided at the lower end by a simple 
drag cable support. The natural frequency of the lateral motion was 
controlled by adjusting the length of a spring rod attached at the 

base of the I beam. Figure 3 presents the variation of model resonant 
frequency with spring rod length. Extreme care was used in the design 
of model components in order to assure repeatability of results and the 
avoidance of model slippage. 

For information on the complete design of the aeroelastic model used 
to represent the antenna boom boards, see Jesberg [71. Included in £7l is 
a description of the displacement transducer used to measure lateral 
motion at the model base, the data handling and tape recording systems, 
the method used to convert analog to digital data, and the digital com- 
puter program employed to convert voltage time history into auto-correla - 
tion and power spectral density forms. 

A listing of the run schedule is provided in Table I. Runs 22 
through 45 were part of the first phase of testing and took place during 
the period of 22 October to 2 November 1965. The objectives of the first 
phase may be stated as: 

RUN OBJECTIVE 

20 to 27 Familiarization and check out of test set-up with 
square and circular cylinders. 

28 to 30 Check the effect of rounding the corners on a square. 
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31 to 36 Investigation of miscellaneous configurations. 

37 to 45 Check the effect of rounding the corners on a square. 

The natural frequencies used in phase one were 36, 42, and 82 cycles 
per second (cps) which corresponded to beam lengths of = 6.5, 5.0, 

D 

and 3.0 inches respectively. Upon completion of the first phase of 
testing, the data that had been recorded on magnetic tape were reexamined. 
Frequencies of all indicial response checks and tunnel test data were 
confirmed by plotting a portion of the tape data playback onto an ink 
recording oscillograph . 

Following phase one, modifications were made in order to improve 
data acquisition and reliability. The spring rod support bracket was 
redesigned for ease of s-etting spring rod length and for improvement of 
spring rod constraint. The displacement sensor was relocated to the 
centerline of the aluminum I beam. In addition, the spring rod was case 
hardened . 

The second phase of these investigations began with run 46 on 24 
January and continued through run 84 on 25 February 1966. The objectives 
for this phase of the program were: 



RUN 

46 to 51 



OBJECTIVE 



Check the effects of modifications and compare 
results with first phase data. 

Check shroud for retrofit on currently installed 
boom boards . 

Check scaling laws and effects of H/D ratio. 

Check feasibility of using a shroud on flanges 

of the circular boom board. 

Investigate the forcing function. 

The natural frequencies used during this phase were different due 

to the equipment modifications that took place between the first and 

second phases. The three basic natural frequencies used were 55, 70, and 

89 cps, which corresponded to P values of 6.5, 5.0, and 3.0 inches 

B 

respectively . 



52 to 62 

63 to 74 
75 to 77 

78 to 84 
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3. Scaling. 

Accurate scaling of a wind-tunnel model is a fundamental considera- 
tion in establishing an experimental program. Important variables in- 
volved in an aeroelastic simulation include fluid density, viscosity, 
velocity, model vibration frequency, stiffness, and mass distribution. 
Since the model was restricted to a single degree of freedom (in the 
lateral direction), the moment of inertia did not have to be scaled. 

All essential scaling variables are included in either Reynolds 
number (Re = VH/-/) or Strouhal number (S = fH/V). The Reynolds number 
of the model and prototype were matched to assure that the aerodynamic- 
forcing functions corresponded. All models were operated at values of 
Reynolds number between 250,000 and 750,000. The full-scale antenna 
element with a frontal width of one foot would have a Reynolds number of 
700,000 for a wind velocity condition of 75 miles per hour (mph) at sea 
level. 

Strouhal number is the dimensionless vortex-shedding frequency of 
the air flow and a match of this number assures a correspondence between 
the model and the prototype dynamic response. In general the prototype- 
antenna elements have a natural frequency of about five cycles per second; 
and as suggested by USNCEL personnel, a dynamic matching should correspond 
to a wind velocity of about 75 mph on a beam with approximately one-foot 
frontal width. At the design point, the fundamental structural Strouhal 
number is about 0.05. Since the model indicated large amplitude responses 
at S = 0.14 and multiples above, it is unlikely that the boom-board 
natural frequencies could act to reinforce the amplitude response at wind 
velocities above about 30 mph. 

It is customary to present wind tunnel results in the form of dimen- 
sionless coefficients, i.e., lift coefficient (C^), drag coefficient 

Li 

(Cj^), etc. In contrast the dynamic responses reported here are not 
normalized and in general it will be necessary to make relative compari- 
sons to see improvements or changes due to configuration revisions. 

It is interesting to check for dynamic matching since a good indi- 
cation of dynamic model characteristics can be obtained. Figure 4 has 
been arranged to show the lateral response of the model for the case of a 
square cross section. Parameters shown in Figure 4 include model size 
(H) and spring rod length (Q^). A casual glance shows that shortening 

D 
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the spring rod reduces the dynamic response, but decreasing the frontal 
width does not consistently decrease the motion amplitude. Therefore, 
comparisons can be made only when taking into account the scaling 
parameters . 

It can be shown that the root-mean-square (R,M,S.) of the lateral 
response is proportional to the expression below for a single degree of 
freedom model : 



\fr ^ fcT ^ 



HjJ 

X 



d>L (f^ 








where : 

y - lateral response at the lower end of the model 
H - beam frontal width 
L - beam exposed frontal length 

I - system mass moment of inertia about the flexure pivot 
f - frequency, cps 

- mean square of the fluctuating lift, assumed to be uniform over 
model length 

spectral density of input aerodynamic forcing function 
Now by a basic definition. 



CJ = 1, cj>i C4") di = £ tjji clS 



( 2 ) 



Equation 2 provides a scaling relation when (|) is expressed in terms of 

density with respect to either frequency or Strouhal number. Since 

f/f = S/S at a given operating condition, the scaling law is: 
n n 
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This law holds providing that and C . - C 

If the Reynolds number effect can be neglected and if there were no 
aeroelastic feedback dependence, then one could assume 'Tc!? -nTcC 
On this basis equation 3 will simplify further to: 



(3) 



n[ _ 9ri Hi / fn 

Njlf 




16 



Comparison again will be made using the data on Figure 4 under the 
assumption of matching Strouhal numbers. The steady-state drag data as 
shown on Figure 8 lead^. to the assumption that Reynolds number effects 
are relatively unimportant for a sharp-edged square or rectangular cross 
section. The starting point for the comparison is from run 58 at q = 20 
psf. The reference values are listed in Table II, part A. In making a 
ratio comparison, the response with the smallest displacement value 
(case c) will be used as a reference. On this basis, calculations were 
made and are listed in Table II, part B. It will be noted that the 
ratios of R.M.S. displacements do not agree with the ratios based upon 
q, H, and f^, but in all cases the displacement ratios were larger. A 
glance at equations 3 and 4 immediately points out that the assumption of 
constant \JcV assumed in equation 4 is probably not valid. Column 3 
in Table II, part B is the ratio showing the discrepancy and has the 
physical interpretation of being the increase in \fc^ due to aeroelastic 
feedback or dynamic motion. 

Figure 5 is a plot showing the apparent increase in fluctuating lift 
due to motion. The abscissa represents lateral motion normalized with 
respect to frontal width. The data points are the values obtained from 
Table II, part B and all values are in reference to an elastic condition. 
A simple rotation of the curve, so that with zero motion the ordinate 
has a value of unity, provides an estimate of the lift increase due to 
aeroelastic feedback. It should be noted that in estimating this curve, 
no information has been obtained as to the magnitude of the fluctuating 
lift. 

A later section of this paper will deal with estimating the 
fluctuating lift and an account will be given to aeroelastic feedback 
again. In the discussion to follow concerning cross sections applicable 
to the antenna problem, only relative comparisons will be made. 
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4. 



Antenna Applications. 



H/D Ratio 

All but one of the boom-board cross sections for the Wullenweber 
antennas are rectangular in shape. The H/D ratio varies from one to 
two and a half. 



WIND 
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A series of runs were conducted in order to demonstrate the influ- 
ence of varying D while holding H constant. The three beams selected for 
this evaluation were: 

(1) 6.0 X 3.0-inch rectangle (H/D = 2.0) 

(2) 6.0 X 4.0-inch rectangle (H/D = 1.5) 

(3) 5.7 X 5. 7 -inch square (H/D = 1.0) 

The small error introduced when using the 5.7-inch square was considered 

acceptable since the values used for comparison were corrected to similar 
values of Strouhal number. 

The responses for natural frequencies of 89 and 70 cycles per second 
are shown on Figures 6 and 7. Table III is the numerical results of the 
H/D ratio study. It will be noted that the amplitude response is approxi- 
mately inversely proportional to the H/D ratio (holding H constant). 

The physical meaning of these results is that the fluctuating 
pressures acting over the lateral faces of the bluff body are primarily 
controlled by the corners and not the width of the lateral faces. 

Corner Radii 

A frequent statement made for circular cylinders is that the unsteady 
lateral air load is of the same order of magnitude as the steady state 
drag while the unsteady drag is an order of magnitude less. Using this 
as a clue, insight into the behavior of circular, square, and rectangular 
cross sections may be gained by looking at Figure 8 which is a summary of 
a drag study made by Delany and Sorensen [3 ] . The square cylinder does 
not show a drag transition trait as for a circular cylinder for Reynolds 
number values well beyond 1x10^. However, placing a radius on the corner 
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edges does bring a transition trait into being. Since the circular-cross - 
sectioned aeroelastic model was relatively more docile thah the square 
cross -section , 5 the material of Delany and Sorenson provided the motiva- 
tion to try corner radii. 

Figure 9 shows a series of cross sections that gradually transform 
the square cylinder to a circular cylinder. The corners of the 5. 7 -inch 
square were successively rounded in one-half -inch increments of radius up 
to and including a two-inch radius. This series of cross sections trans- 
formed the square cylinder into a circular cylinder in definable steps. 

It was anticipated that rounding the corners of a square would produce 
responses that fall between the responses of the square and circular 
cylinders . 

Typical results obtained are the responses for a model with a three 
inch spring rod length as shown in Figure 10. The curve farthest to the 
left, is the response for a square and it will be noted that the response 
rapidly exceeds the limit of the graph. The addition of one-half -inch 

radius corners (R = 0.5 inch) on the 5.7-inch square (H = 5 . 7 inch) 
results in some improvement (R/H = 0.0877). The one-inch radius 
(R/H = 0.176) shows marked improvement. The one and a half -inch radius 
(R/H = 0.263) exhibited a low amplitude response characteristic and the 
two-inch radius (R/H = 0.351) was comparable to the circular cylinder. 

The advantages that can be obtained by rounding the corners of the 
boom boards are obvious when looking at Figure 10. The solid laminated 
boom boards would be particularly well adapted to this fix. However the 
boom boards with hollow centers would be seriously weakened if the corners 
were rounded. 

R/H = 0.263 is suggested for the construction of new boom boards. 
However, a corner radius ratio of R/H = 0.176 would provide a considerable 
improvement over the square corners which are now being used. 

Gloves 

Rounding the corners of a rectangular beam provided well defined 
decreases in lateral amplitude response of the beam. There are two 
reasons that discourage the use of this fix on the currently installed 
boom boards. First, all boom boards that are hollow in the center would 
be seriously weakened if the corners were rounded. CSecondly, the boom 
boards are physically located 100 feet above the ground. It would be 
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easier to install some additional piece of material to a beam than remove 
part of the beam, 

A shroud on the corners of a square was selected as a possible fix. 
Figure 11 is a sketch of this shroud or glove as it will be referred to 
in this paper. From the results of an earlier series of tunnel runs 
where the corners of a square were rounded in successive one -half -inch 
increments, a ratio R/H = 0,22 was selected for the glove. The square 
with a two-inch radius (R/H = 0,351) had the most effect in reducing the 
dynamic response when compared with the basic square cross section. 
However, a glove with an R/H ratio of 0,351 would be relatively large, 
thus greatly increasing the weight and frontal area of the beam, A 
compromise was made in selecting a glove radius so that undue penalties 
in frontal area and weight addition would not take place. wKSildc still 
providing the benefits of reducing dynamic response of the boom boards. 

Two sets of gloves (R = 0.875 inches and R = 1.25 inches) were 
constructed to be used with two squares (H = 4.0 inches and H = 5.7 
inches). This provided four possible configurations with gloves to check 
for possible scaling laws and to give some indication of a desirable R/H 
ratio . 

(1) 5.7-inch square with l,25-^‘,inhh gloves (R/H = 0.22) 

(2) 5. 7- inch square with 0.875-inch gloves (R/H = 0.15) 

(3) 4.0-inch square with 1.25 - inch gloves (R/H = 0.31) 

(4) 4.0-inch square with 0.875-inch gloves (R/H = 0.22) 

Configuration (2) would indicate if the R/H ratio selected were too large 
Configuration (3) would likewise indicate if the R/H ratio were too small 
All four configurations have been plotted in Figure 12 for a spring rod 
length of six and one half inches (f^ = 55 cps ) . All curves exhibited 
considerable improvement over the response of the square cross section. 

Figures 13 and 14 show the response of the 5. 7 -inch squares with and 
without gloves. A ratio of R/H = 0.22 has a clear advantage of a lower 
response below a dynamic pressure of 30 pounds per square foot. When a 
natural frequency of 55 cps was used the R/H =0,15 demonstrated a larger 
peaked displacement at 10 pounds per square foot dynamic pressure. 

Figures 15 and 16 show the response of the four-inch square with and 
without gloves. R/H = 0.31 and R/H = 0.22 are relatively the same 
amplitude at f^ = 70 cps (Q^ = 5.0 inches). However, at f^ = 55 cps 
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(J^ =6.5 inches) R/H = 0.22 has a definite advantage below a dynamic 

6 

pressure of 30 pounds per square foot. 

It is clear that the best radius for rounding the corners of a 
square is definitely not the best radius to be used for gloves on that 
same square. The ratio R/H when using gloves should be about 0.22. 

Circular Cylinder 

A positive advantage was derived by shrouding the sharp corners of a 
square. Therefore a brief look was taken at shrouding spanwise flanges 
located at the maximum width points (©=^9(5) on the circular cylinder. 

The recommended fix for the circular cylinder by Hite and Breckon 
L2] was the rotation of the beam so that the flanges would be located at 
the stream stagnation points. Utilization of this fix in the field was 
infeasible for various physical reasons. 

Hite and Breckon demonstrated that the flanges when located at 
were responsible for the large oscillations of the circular cylinder. 
Without the flanges the flow separated randomly aft of the 90 degree 
position (present location of the flanges on the circular boom board). 

With the flanges installed, separation of the flow was stabilized to occur 
at the flanges, concentrating the energy of the lateral force at a 
dominant Strouhal number. The concentration of random energy in the 
frequency spectra results in a violent response. Therefore if the 
flanges can not be removed or rotated to stream stagnation points, the 
sharp edge effect of the flange should be eliminated. 

Shrouding of the boom-board flange was expected to produce desir- 
able responses. A one-inch radius shroud was selected, which was a ratio 
of shroud radius to beam diameter of 0.198 (R/H = 0.198). The shroud was 
extended laterally (into the stream flow) only far enough to cover the 
flange. This was a distance of only one quarter of an inch on the model, 
thus the cylinder retained a relatively smooth shape as is illustrated in 
Figure 17. The shroud was glued to the front half of the beam to provide 
the smoothest possible surface for the air flow. 

The responses obtained are plotted on Figures 18 and 19. With a 
spring lever arm of three inches, the response of the possible fix 
exceeded the response of a square (approximately the result of a flange 
at 90 degrees) at a dynamic pressure of 15 pounds per square foot. The 
proposed fix will apparently also exceed the square's response at about 
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10 pounds per square foot when a spring rod length of six and one-half 
inches is used. A larger shroud radius relative to beam diameter might 
give better results. Another solution might be to blend the shroud 
fairing smoothly into the cylinder contour without discontinuity in 
surface slope. 

Further expenditure of effort was not considered justified because 
the circular boom boards at Winter Harbor, Maine have suffered further 
deterioration. Replacement of the circular cylinder with rectangular- 
cross -section boom boards is under consideration. 
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5 . 



Estimate of Unsteady Lift. 

The last six runs (78 to 84) of the program were conducted in order 
to gain some knowledge of the lateral forcing function. Both the co- 
efficient of lift and its spectral distribution were estimated. A sharp 
edged configuration was avoided for this evaluation so that a full set of 
data could be obtained without restriction on amplitude response. For 
this reason, the 5.7-inch square with one-inch corner radii (R/H = 0.176) 
was selected. The natural frequency of the beam was varied from 60 to 
89 cps by adjusting the spring rod length (J?„) in one-half -inch incre- 

D 

ments from six to three inches. Data were obtained at tunnel dynamic 
pressures up to 40 psf using intervals of two psf. 

The amplitude responses (recorded on magnetic tape) were digitized 
to permit the use of a high-speed digital computer [7 ] . Utilizing the 
auto-correlation function, the power spectral density (P.S.D.) of the 
system output response was computed by the Tukey -Blackman method [l] . 
Sample power spectral densities of the output are shown on Figures 20 
and 21. These curves are computer plotted and represent the complete 
set of values at dynamic pressures of eight and 40 psf respectively. 

It will be noted on Figure 20 that there are energy peaks at 24, 50, and 
80 cps. These correspond to Strouhal numbers of 0.14, 0.28, and 0.45 
respectively. The large energy peak at S = 0.14 was evident for values of 
q up to about 15 psf (Re = 340,000). The characteristics shown on 
Figure 21 at q = 40 psf (Re = 560,000) indicate energy peaks at f = 41, 

80, and 123 cps which correspond to Strouhal numbers of 0.11, 0.21, and 
0.32 respectively, and in general the majority of energy is at the higher 
Strouhal numbers. These traits were shared by the major portion of out- 
put response data when Reynolds number was above 340,000. 

Computation of the power spectral density (P.S.D.) of the output 
from recorded data opens the door to investigation of the forcing 
function. 4^o(f^ the P.S.D. of the output response was normalized. 



The upper limit of 125 cycles per second assumes that there is no signi- 



this assumption was checked and found to be correct for the calculations 
used in this report [7] . 




(s) 



ficant contribution to 




The accuracy of 



is the P«S.D, of the input forcing function, (a distributed 
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force acting along the length of the model, but which may be treated as 
a concentrated load acting at the mid-span of the model). is 

related to through the following relation. 

X (;4') is the transfer function which is dependent only on the model 
being tested. 

l2a<^")l^ = ([i-o-(x) ("■!■?] (j) 

f is the frequency in cycles per second and f^ is the natural 
frequency of the beam. ^ is the damping ratio of the beam. It was 
calculated by log decrement procedures on indicial-response checks. Due 
to the difficulty of reading the amplitudes of successive peaks with 
sufficient accuracy, an average ^ of 0.018 was used for computation of 

. This value was obtained from eleven log decrement calculations. 
The damping ratio was sensitive to spring rod clamping. Various values 
of ^ (0.018 to 0.027) were used in calculating (pX (4) to check the 

effect of ^ • A 50 per cent increase in ^ resulted in less than one 

per cent change in 

The forcing function may be found by considering the geometry of 
the test set up in the wind tunnel as sketched below: 




Values used were: 

Xcc = ^-0 inches 

X =46.0 inches 
Xc =57.0 inches 
Xj =54.0 inches 
I = 3.5.73. slug- in^ 



The equation of motion for harmonic excitation may be expressed as: 



le + ce -t- KXc© = F s\n tot 



( 8 ) 
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Replacing the angular coordinate by the displacement at the transducer 
yields : 



Y * Xa 6 

s,n U.t 

Y-^Ei^u^Y + t-o^Y = Xb + x<x j 

let M - I L 2 J 

/s 

V = , ■ 

— ^ 

If the mean square of the displacement is Y j then the mean square of 
the applied load is: 



where 



J<L f Xb-^-XcL ~|^ 

11 = ja L E. J 

( M-'h') ~ Jo (f)o(-f) 



= (M.s.Uy^) 

I 2 I cL-f 



and 



iF = 



±d. 







H CXb-Xa.) 



therefore 




(ex) 

<3 H Xa ( Xt -XI) 






S/Ct is the root -mean-square of the coefficient of lift. 

A high-speed digital computer was utilized to calculate the root- 
mean-square coefficient of lift for each of the seven natural frequen- 
cies of the beam and at each dynamic pressure used. Figure 22 is the 
plot of the individual >TcF and Figure 23 is a plot of the average 
Both plots clearly indicate the occurrence of a transition at 14-pounds 
per square foot dynamic pressure (approximate Reynolds number of 
3.3 X 10^) where the flow exhibits a change to supercritical flow. 

All curves are essentially parallel and with only minor exceptions the 
beams with lower natural frequencies have larger coefficients of lift. 
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The consistent pattern of these plots suggests that aerodynamic forces 
may be induced by the motion of the models through an aeroelas tic-feed- 
back influence. 

The P.S.D. of the input was also machine plotted using the digital 
computer to handle the transfer function type calculations. Figures 24 
and 25 represent estimates for subcritical and supercritical flow and 
show good agreement for all values of model resonant frequency '^rith 
respect to the occurrence of energy peaks. 

The peaks of the power spectral density curves of the exciting 
force occurred at a constant Strouhal number. For dynamic pressure up 
to and including 12 pounds per square foot (Re = 3.1x10^) the vortex 
shedding occurred at a Strouhal frequency of 0.140. At and above 16 
pounds per square foot (Re = 3.5x10^) the Strouhal frequency was 0.32. 

At 14 pounds per square foot (Re = 3.3x10^) two peaks were present, 
one at a Strouhal number of 0.14 and the other at 0.32. The beams with 
natural frequencies below 75 cycles pertsecond (70 65, and 60 cycles 

per second) responded at a Strouhal number of 0.32. 

The 14 pounds per square foot reading was obviously in the flow 
transition region. Supercritical flow was experienced at a dynamic 
pressure of 16 pounds per square foot (Re = 3.5x10^). Dr. Fung found 
that the critical Re 3 molds number for his circular cylinder was 
Re = 3 .0x10^ [6 ] . 

The aeroelastic— feedback term was further investigated by consid- 
ering results of Figure 22 at dynamic pressures of 10 and 30 psf 

which correspond to Reynolds numbers of 280,000 and 480,000 respectively. 
The fluctuating lift values are cross plotted on Figure 26 as a function 
of the normalized lateral response. It will be noted that both flow condi- 
tions have approximately linear dependence upon motion, and the sub- 
critical flow appears to extrapolate to a value of nTcT = 0*35 for a 
rigid model while the corresponding supercritical flow value is about 
0.075. A graphical procedure has been used on Figure 26 to obtain an 
aeroelastic gain correction for the subcritical flow case, normalized 
such that the gain value is zero for the rigid case. The large influence 
of motion on fluctuating loads is quite surprising, especially when one 
notes that lateral displacement (R.M.S.) on the order of one hair thick- 
ness more than doubles the value of fluctuating load. 
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A simple model may be conjectured to describe the build-up of an 
aeroelastic feedback term. The model assumes that fluctuating lift is 
developed due to the instantaneous response of a rotating drag vector, 
with the rotation being caused by the motion-induced angle of attack 
changes. Consider a bluff body moving laterally in rectilinear flow as 



shown in 


the sketch below: 








r 


If 


y = y sinwt, 


then 


= y = Vo sin 


hence 


— S.'TT ‘f’ Vo 


and 





Therefore the mean-induced angle is given by: 






= 






V‘ 



But if there exists a dominant Strouhal number that defines the frequency, 
then f = VS/H. 



Hence = ETT®" 

If we assume that 



then 



and 



Cu ('fc') — ^ (‘fc') C-O steady 
D steady 



nTF =vr?iTS^c, 






steady 



(< 0 ) 

( 11 ) 



For the square cylinder with rounded corners under the condition of 
subcritical flow, we may assume S = 0.14 and = 1.0. Hence an estimate 
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of the aeroelastic effect would be 



= nTF 'TT (o.l4)(|.o') = 6.2. 10*' 

® vTh ) 



The actual slope estimated from Figure 26 is 5x10 , a value 
approximately four orders of magnitude greater than that predicted by 
the simple model. The large discrepancy would indicate that either the 
experimental procedures and estimates presented here are suspect of that 
the simple model is completely inadequate. 

The existence of fluctuating loads on a rigid model seems reasonable 
since the unsteady pressure exist in separated flow, independent of motion. 
The dependence of these pressures and their instantaneous integrated 
value (net fluctuating force) upon motion is not clear, and further study 
certainly seems worthwhile. 
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6. Summary. 

The large responses of the square-cross -section boom boards were 
sharply reduced by two different modifications to the beams. Rounding 
the corners on solid boom boards at the time of construction appeared to 
be a desirable solution. A ratio of corner radius to frontal width of 
0.263 provides a response that is near the level of a circular cylinder 
for all wind velocities (up to 171 miles per hour) that may be encountered. 

A retrofit modification was tested for use on boom boards that are 
presently installed as part of the Wullenwebber Antenna system. This fix 
consisted of shaped shrouds which would be attached over the corners of 
the boom boards. The desirable ratio of shroud radius to the length of 
the side of the square should be ^boutC0.22Lc C.r. . 

A shroud for the flange on the circular cylinder was tested, but no 
appreciable improvement was noted. A different size and or shape of 
shroud might produce a marked reduction in amplitude response, but fur- 
ther studies in this direction were not made. 

Tunnel tests demonstrated a correlation between models of different 
H/D ratios (ratio of frontal width to depth of model) in the Reynolds 
number supercritical flow region. The first order components of the 
amplitude responses were found to be inversely proportional to the H/D 
ratios (when holding frontal width constant), thereby showing that the 
fluctuating pressures developed on the lateral faces are not too sensi- 
tive to the length of that face. The dominant factor appears to be 
corner radius ratio. 

The use of precision instrumentation, magnetic recording of data, 
and data processing by analog and digital computers permitted investi- 
gation of the forcing function's power spectral density and lateral 
root -mean -square coefficient of lift. The coefficients of lift were 
estimated for the square with rounded corners and were compared with 
previously published information on circular cylinders. The dependence 
of unsteady lift upon motion strongly suggests the influence of an 
aeroelastic feedback term. Further investigations should be conducted 
in this area with different size models to confirm and/or isolate the 
aeroelastic feedback component. 
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35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 
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RUN SCHEDULE 



Configuration 
Circ Cyl 
Circ Cyl 
Circ Cyl 
Circ Cyl 
Circ Cyl 
Sqr Cyl 
Sqr Cyl 
Sqr Cyl 

5.7 Sqr & .50 RC 

Front-5.7 Sqr & .50 RC 
Rear-Rect Cyl 

I 

Rect Cyl 

I 

Front -Rect Cyl 
Rear-5.7 Sqr & .50 RC 

Sqr Cyl & 1.00 RC 

V 

5.7 Sqr & 1.50 RC 

5.7 Sqr & 2.00 RC 

'r 



q(PSF) 
20 
5-10 
5-40 
5-40 
5-40 
5-20 
5-12 
3- 8 
3-10 
5-15 
5-30 
5-25 

3-10 
3- 8 
5-15 
5-20 
5-15 

5-40 

u 

5-40 



5.7 Sqr & 2.00 RC 
5.7 Sqr & 1.25 Glove 



5-:40 



8.88 

8.88 

8.88 

6.00 

3.00 

3.00 

5.00 
6.50 
6.50 

5.00 

3.00 

3.00 

6.50 

6.50 

5.00 

3.00 
3.00 

3.00 

5.00 
6.50 
6.50 

5.00 

3.00 

3.00 
6.50 

5.00 
5.00 
3 .00 
5.00 
6.50 



Data 

F,A,D 

F,A,D,T 



D 

D 

D 

D 
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Run 

sF" 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 
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RUN SCHEDULE (CONT'D.) 



Configuration 



q(PSF) X 



E_ 



5.7 Sqr & 1.25 Glove 
5.7 Sqr & 1.25 Glove 

T 

5.7 Sqr & .875 Glove 



5,40 



4.00 Sqr 1.25 Glove 



I 



6.00 X 4.00 Rect 



6.00 X 3.00 Rect 



Giro Cyl 1.00 R over 
flange 



5.7 Sqr & 1.00 RC 

I 



3- 8 
3-15 
5r25 
5.25 
5r20 
3 t10 
5-40 

5t40 

5-35 

5.40 

t 



6.50 

3.00 

3.00 

5.00 
6.50 
6.50 

5.00 

3.00 



5.7 Sqr 


5-20 


3.00 


5.7 Sqr 


5-15 


5.00 


4.00 Sqr 


5-12 


5.00 


4.00 Sqr 


3- 8 


6.50 


4.00 Sqr 


5-20 


3.00 


4.00 Sqr & .875 


5-40 


3.00 



5.00 

6.50 

5.00 

3.00 
6.50 
6.50 

5.00 

3.00 

3.00 

5.00 
6.50 

6.50 

5.00 

3.00 
3.00 

3.50 



Data 

A,D 

A,D 

F,A,D,T 



r 

F,A,D,T 



F,D,T 
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RUN SCHEDULE (CONT'D.) 



Date Run Configuration 

2-25 80 5.7 Sqr & 1.00 RC 

81 
82 

83 

84 r 



q(PSF) 

5-40 



t 



4.00 

4.50 

5.00 

5.50 

6.00 



RC - Radius of Corners 
Circ - Circular 
Cyl - Cylinder 
Sqr - Square 
Rect - Rectangular 

- Length of Spring Rod 



Data Code: 

F - Frequency Response Measured (Tunnel-off) 
A - Acceleration Readings Recorded 
D - Displacement Readings Recorded 
T - Data Recorded on Magnetic Tape 



All lengths in inches 



Data 

F,D,T 



ir 
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Table II 



Scaling Study Calculations 
Square Cross Section 



A. Reference Data 



Case 


Run 


H 

(in.) 


(in. ) 


f 

(cps) 


q 

(psf ) 


S 

n 


(lO’^in. ) 


(xlO^^) 


a 


58 


5.7 


3.0 


89 


20.0 


0.32 


6.57 


1.15 


b 


59 


5.7 


5.0 


70 


13.0 


0.32 


11.15 


1.96 


c 


62 


4.0 


3.0 


89 


10.0 


0.32 


1.37 


0.34 


d 


60 


4.0 


5.0 


70 


6.5 


0.32 


2.75 


0.69 


e 


61 


4.0 


6 .5 


55 


•',4 1.0 


0.32 


3.20 


0.80 



B . Comparison 



Cases 




f 3. H. irin.f 
© 


(D/© 


(d):(c) 


2.00 


1.05 


1.90 


(e):(c) 


2.34 


1.05 


2.35 


(a):(c) 


4.80 


2.83 , 


1.69 


(b):(c) 


8.15 


2 .99 


2.72 
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Table III 



H/D Calculations 



f 

(cps) 


q 

(psf ) 


H/D 


{T 

(10 ^ inch) 


Amp* 

Ratio 


Inverse 

H/D 

Ratio 






2.0 


4.79 


.70 


.75 




25 


1.5 


6.85 


1.00 


1.00 






1.0 






1.50 


89 




2.0 


3.21 


.67 


.75 




20 


1.5 


4.79 


1.00 


1.00 






1.0 


6.57 


1.65 


1.50 






2.0 


1.97 


.73 


.75 




15 


1.5 


2.70 


1.00 


1.00 






1.0 


3.90 


1.44 


1.50 






2.0 


9.43 


.67 


.75 




15 


1.5 


14.0 


1.00 


1.00 






1.0 


20.2 


1.44 


1.50 






2.0 


4.79 


.75 


.75 


70 


10 


1.5 


6.4 (EST) 


1.00 


1.00 






1.0 


8.6 (EST) 


1.35 


1.50 






2.0 


2.03 


.77 


.75 




5 


1.5 


2.63 


1.00 


1.00 






1.0 


3.60 


1.37 


1.50 



* Amplitude ratio is defined relative to the amplitude of the bluff 
body with H/D = 1.50. 
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FIGURE 1 

ANTENNA INSTALLATION 

general layout 
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FIGURE 2 

TYPICAL BOOM BOARD 
CROSS SECTIONS 








7 


z 







- 8 "-] 



— 8 *- 
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FIGURE 3 

SPRING BAR CALIBRATION 
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FIGURE 4 

DYNAMIC RESPONSE, SQUARE CYUNDER 




I 
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FIGURE 5 

EFFECTS OF MOTION ON LIFT, 
SQUARE CYLINDER 



Data From 
FUN to RUN 
0 62 62 

A 60 * RUN 62 Is Elastic Ref. Value. 

S'! " See TABLED for Computation. 

V 58 • 
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FIGURE 6 

EFFECTS OF H'/D RATIO, ig=3 INCHES 



RECTANGULAR CROSS SECTN. CYLS., J! =3 OIn. 

D 



SYM 

A 

□ 

© 



D ^ 






1.0 58 y 






1.5 71 




Hs 

1 


2.0 72 







:6 in. 




Dynamic Pressur* (#/ft2) 
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Lateral 
Displ. 
-RMS- 
(10*3 In.) 



FIGURE 7 

EFFECTS OF H/t) RATIO, 



S„-5 INCHES 

o 



RECTANGULAR CROSS-SECT. CYLS., 



SYM 


H 

D 


RUN 


A 


1.0 


59 


□ 


1 .5 


70 


O 


2.0 


73 




J?g=5.0in. 









4 

H = 6 in 

J. 
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FIGURE 8 

BLUFF BODY DRAG SUMMARY 




OVUO JO 1N3DIJJ30D 



43 



REYNOLDS NUMBER (10^) 



FIGURE 9 

MODEL CROSS SECTIONS 



^5.7" - 






t 

5.7 

1 



Square 



r 







-j 



0.5" Radius 
R/H =0.0877 






^ 

1.0” Radius 
R/H =0.176 



\ 



7 

1.5" Radius 
R/H=0.263 




J 

2.0" Radius 
R/H= 0.351 




Cylinder 
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FIGURE 10 

EFFECTS OF CORNER RADII 
ON DYNAMIC RESPONSE 



o 
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DYNAMIC PRESSURE (*/FT*) 




FIGURE II 



GLOVE GEOMETRY 



WIND 
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FIGURE 12 
EFFECTS OF GLOVES 
^8 »6.5 INCHES 




(HDNit.oixi) aonindwv 



47 



DYNAMIC PRESSURE (*/FT®) 



FIGURE 13 

EFFECTS OF GLOVES 
ON 5.7 -INCH SQUARE 
Pb ’6.5 INCHES 



# 




(honi j.oixi) aofundwv 
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DYNAMIC PRESSUF^E (#/FT*-) 



FIGURE 14 
EFFECTS OF GLOVE 
ON 5.7 -INCH SQUARE 
Pa = 5.0 INCHES 




«0 ^ 



CH3NI jjoi i) aonindwv 
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DYNAMIC PRESSURE (#/FT® ) 



FIGXmE 15 

EFFECTS OF GLOVES 
ON 4.0 -INCH SQUARE 
Jb= 6.5 INCHES 




(HDNi s.oixi) sonindwv 
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DYNAMIC PRESSURE (^f/FT*) 



FIGURE 16 
EFFECTS OF GLOVE 
ON 4.0 -INCH SQUARE 
Jb = 5.0 INCHES 




(honi s.oixi) aanindwv 
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DYNAMIC PRESSURE (# /FT^) 



FIGURE 17 



CIRCULAR CROSS SECTION , 
SHROUD FAIRING 
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FIGURE 18 

EFFECTS OF SHROUD FAIRINGS 
CIRCULAR CYLINDER ^b"30 INCHES 




(.HDNi ,.oixi) 3arundwv 
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DYNAMIC PRESSURE (#/ft‘) 



FIGURE 19 

EFFECTS OF SHROUD FAIRINGS 
CIRCULAR CYLINDER Pe=6-5 INCHES 




(HDNI j.O|x I) 3onindwv 
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DYNAMIC PRESSURE (*/fT*) 
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Figure 20 

P.S.D. OF OUTPUT AT q=8 PSF 
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Notes : 

1. Square cylinder and corner 
radii, R/H=0.176 
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9 
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Figure 21 

P.S.D. OF OUTPUT AT q=40 PSF 

I 




Frequency (cycles/second) 
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REYNOLDS NUMBER (10^) 



FIGURE 23 
COEFFICIENT OF LIFT 
DATA SUMMARY 
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REYNOLDS NUMBER 



010 015 020 025 ^ 



Figure 24 

P.S.D. OF INPUT AT q=8 PSF 







Notes: 

!• Square cylinder and corner 
radii, R/H= 0.176 
2. Re= 5.6 X 10^ 
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Figure 25 

P.S.D. OF INPUT AT q=40 PSF 



Notes: 

1. Square cylinder and corner 
radii, R/H=0.176 
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FIGURE 26 

AEROELASTIC EFFECTS ON UNSTEADY LIFT 



SQUARE CYL. ♦CORNER RADII; ^*0.176 
SYM ^ ^ 

O 0.28x10®, 10 PSF 
A 0.48x10®, 30 •< 
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been modified such that the energy spectra has apparently been broadened from a 
single dominant peak into many weaker peaks. Thus the effective total response 
is reduced. 

The digital data processing has been extended in scope due to the repeat- 
ability in signal analysis. Estimates have been made of the "two dimensional" 
random input forcing function using the output response in conjunction with the 
aeroelastic model transfer function. Reasonably good results were obtained and 
an estimate was made of the root mean square value of the unsteady fluctuating 
air load. 
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